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SUMMARY 


"~-^A  solution  is  presented  for  the  elastic  stress  intensity  factors  at  the 
tips  of  a  slightly  kinked,  partially  closed  crack  in  compression.  The  solu¬ 
tion  is  accurate  to  first  order  in  the  deviation  of  the  crack  surface  from  a 
straight  line  and  is  carried  out  using  perturbation  procedures  analogous  to 


those  of  Banichuk  [5] [  Goldstein  and  Salganik‘[6]  and  Cotterell  and  Rice  T7j 

)  V 

for  the  problem  of  an  open  crack.  Comparison  with  the  exact  solution  indi¬ 
cates  that  the  asymptotic  solution  is  accurate  for  values  of  the  angle  between 


the  straight  crack  and  its  out-of-plane  kinks  up  to  about  20  . 


ir 


1.  INTRODUCTION 


Experiments  on  glass  plates  containing  pre-existing  planar  through  cracks 
oriented  at  an  angle  to  the  direction  of  the  axial  compression  have  revealed 
that  the  relative  sliding  of  the  faces  of  the  pre-existing  cracks  does  not 
result  in  co-planar  crack  growth,  bu'  rather  produces  at  the  tips  of  the  pre¬ 
existing  cracks  small  tension  cracks  which  deviate  at  sharp  angles  from  the 
sliding  plane  [1-4].  These  experiments  are  designed  to  be  models  for  the  pro¬ 
pagation  of  cracks  in  rocks  in  compression.  In  this  paper,  we  are  concerned 
with  the  calculation  of  stress  intensity  factors  at  the  tips  of  the  kinked 
open  extensions  of  a  closed  sliding  through  crack.  The  same  method  can  be 
extended  to  a  curved  crack  with  several  closed  sections.  The  solution  ob¬ 
tained  is  accurate  to  first  order  in  the  deviation  of  the  crack  surface  from  a 
straight  line  drawn  between  the  kink  tips  and  is  carried  out  using  perturba¬ 
tion  procedures  similar  to  those  used  in  Refs.  [5-9]  for  the  problem  of  the 
open  crack.  The  results  can  be  stated  in  terms  of  known  solutions  for  a 
single  straight  crack  or  a  co-linear  array  of  straight  cracks. 


2 


A  complete  solution  to  the  problem  of  the  sliding  kinked  crack  has  been 
given  by  Nemat-Nasser  and  Horii  [3],  who  used  a  continuous  distribution  of 
dislocations  to  model  the  crack  and  its  kinks.  Zn  order  to  find  the  stress 
intensity  factors,  they  solved  numerically  a  singular  integral  equation  for 
the  dislocation  distribution.  In  contrast,  we  can  avoid  the  solution  of  the 
singular  integral  equation  by  using  the  results  of  the  asymptotic  analysis  for 
the  stress  intensity  factors.  However,  the  validity  of  the  asymptotic  solu¬ 
tion  is  limited  to  small  deviations  of  the  crack  surface  from  a  straight  line. 
Comparisons  with  the  exact  solution  given  in  Ref.  [3]  Indicate  that  the  first 
order  solution  for  the  mode  I  stress  intensity  factor  is  accurate  for  values 
of  the  angle  between  the  straight  crack  and  its  out-of-plane  kinks  up  to  about 
20°. 


2.  GENERAL  FORMULATION  OF  THE  PROBLEM 


2.1  Formulation  of  the  boundary  value  problem 

Consider  an  infinite  plate  of  a  homogeneous,  isotropic,  linearly  elastic, 
brittle  solid  containing  a  curved  crack  on  y  =■  \(x)  ,  with  its  tips  at  the 
positions  x  =  ±a  (Fig.  1).  A  uniform  state  of  stress  ,  a ^  and  ax” 

is  applied  at  infinity,  with  Oyy  <  0  and  oXy  <  0  ,  where  tension  is 
regarded  as  positive.  The  corresponding  two-dimensional  boundary  value 
problem  is  given  by 


°jt,j  -  0 

2£ij  "  ui,j  +  uj,i  ; 

I 

; 

I 


°ij  ”  cijklek! 
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onn(x,  X)  -  <Jn8(x,  X)  =■  0 

<Jn8<x»  x>  “  ^nn<x»  x> 
u£(x,  X)  =*  u"(x,  X) 


at  infinity,  (2) 
on  the  open  portions  of  the  crack,  (3) 
on  the  sliding  portions 

of  the  crack,  (4) 


where  Oj j  ,  e^j  and  are  the  stress,  strain  and  displacement  fields  in 

the  region  V  occupied  by  the  body,  Is  the  fourth  order  tensor  of  the 

elastic  moduli,  onn  and  on8  are  the  normal  and  shear  tractions  at  the 
crack  surface,  un  is  the  displacement  in  the  direction  normal  to  the  crack 
surface,  p  Is  the  coefficient  of  friction,  A^  Is  SA/bx^  and  the  super¬ 
scripts  plus  and  minus  denote  the  value  of  the  Indicated  quantity  on  the  upper 
and  lower  surfaces  of  the  crack.  Note  that  the  open  and  sliding  portions  of 
the  crack  are,  in  general,  not  known  In  advance  and  their  determination  be¬ 
comes  part  of  the  solution. 

2.2  Small-parameter  expansion 

The  essence  of  the  approximation  we  use  is  that  the  solution  to  the  pro¬ 
blem  with  the  curved  crack  is  close,  In  some  sense,  to  the  solution  of  a 
similar  problem  for  a  straight  crack.  In  fact,  we  shall  use  the  solution  to 
the  following  problem,  Involving  a  flat  crack,  as  the  leading  or  zeroth  order 
approximation  in  our  expansion.  Let  and  u^^  be  such  that 


,<°) 

ji  >  j 


9  (0)  (0)  (0) 

2e, .  =  u,  ,  +  u,  , 

ij  i.j  j.i- 


ij 


C  c(0) 

ijkJrki 


in  V' 


(5; 


at  infinity 


(6) 


and 


0)  -  «,«>Cx,  0)  -  0 
0)  -  IWyJhx,  0) 

U^X,  0)  “  Uy(X,  0) 


on  the  open  portions  of  the  crack,  (7) 
on  the  sliding  portions 

of  the  crack,  (8) 


where  V'  is  the  plane  with  a  straight  slit  lying  on  the  x-axis  from  -a 
to  a  .  If  the  slope  of  the  actual  c rack,  \'(x)  ,  has  order  of  magnitude 
e  «  1  at  its  largest,  then  we  can  seek  a  perturbation  expansion  in  e  for 
the  solution  to  the  problem  of  the  curved  crack,  such  that 


„  (0)  (1) 

0*0  +  O 

»W  M  /v 

e  -  e<°>  +  e<» 

„ .  +  .U) 

/V  <V  /v 


+  0(e2)  , 
+  0(e2)  , 
+  0(e2)  , 


(9) 

(10) 

(ID 


where  and  ^  are  all  0(e)  compared  to  the  leading  order 

terms.  We  mention  that  X’(jc)  =  0(e)  also  means  that  k(x)/a  «  0(e)  ,  be¬ 
cause  X(±a)  -  0  .  What  remains  now  is  the  finding  of  the  equations  and  the 
boundary  conditions  governing  ,  e^^  and  u^^  .  We  note  at  this 

stage  that  our  approach  is  identical  to  that  of  Cotterell  and  Rice  [7],  except 
that  they  addressed  the  problem  of  a  crack  open  everywhere.  Furthermore,  they 
found  their  solutions  and  expressed  their  expansions  in  terms  of 
Muskhelishvili’s  [10]  complex  potentials.  We  prefer  to  work  in  terms  of  fun¬ 
damental  quantities,  although  it  is  entirely  possible  that  the  partially 


closed,  slightly  curved  crack  can  also  be  solved  by  a  variation  of  the  complex 
variable  treatment  of  Cotterell  and  Rice  [5]. 

We  return  now  to  the  question  of  finding  and  u^^  .  In 

order  to  find  the  equations  and  boundary  conditions  governing  ,  e^ 

and  «<l>  ,  we  substitute  the  expansions  (9)-(ll)  into  equations  (l)-(4). 

We  also  use  the  fact  that  both  X(x)  and  X'(x)  are  0(e)  to  write 
expansions  in  e  for  the  tractions  and  displacements  on  the  crack  surface  y 
■  \(x)  .  Using  a  tensorlal  transformation,  we  find  that  the  normal  and  shear 
tractions  on  the  actual  crack  can  be  written  as 

a  (x,  X)  **  -=•  [<J  (x,  X)  +  a  (x,  X)]  + 

nn  *  2  xx  ’  yy  ’ 

+  J  [<7yy(x,  “  °XX^X’  ^^COS29  “  CTxy^X»  X)stn20  » 

<*ng(x,  x)  “  °xy(x’  X>cos20  +  Y  l®yy Cx*  x)  -  (x,  X)Jsin20  , 

where  9  »  X'(x)  +  0(e)  .  Then,  using  a  Maclaurln  series  expansion  in  0 

for  sin29  and  cos20  ,  we  find 


°nn^x»  x)  “  ayy(x>  x)  **  2X*(x)  axy(x,  X)  +  0(e2)  , 

CTns  (x,  X)  -  axy(x,  X)  +  X'(x)[Oyy(x,  X)  -  axx( x ,  X)]  +  0(e2)  . 


If  we  now  write  Maclaurln  series  expansions  in  y  for  <j  ,  avv  and  avv  , 

y y  xy 

the  last  two  equations  become 


9oxy(x’  0) 


2X r (x)  o  (x ,  0)  +  0(e2)  , 
xy 


a  (x,  X)  ■  a  (x,  0)  -  X(x) 
nn  yy 


(12) 
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a<j  (x,  o) 

^(x,  X)  •  ^(x,  0)  -  X(x)  - ^ - +  X'(x)[cyy(x,  0)  - 

“  <*xx(x,  0)]  +  0(c2)  , 


where  we  have  also  used  the  equilibrium  equations  8oxv/5y  -  -dcrxx/5x  and 


dOyy/dy  ■  -dOx y/dX  . 

In  a  similar  way  we  can  show  that 


u^x,  X)  -  uy(x,  0)  +  X(x)6yy(x,  0)  -  X'(x)ux(x,  0)  +  0(e2) 


Using  the  expansions  (9)-(ll),  equations  ( 12)— ( 14)  can  be  written  as 

/,v  ao(0)(x,  0) 

ffnn<X’  X)  “  °yy)(x’  0)  +  °yy)(x*  0)  "  X(x>  “ - 

-  2X»(x)cr^)(x,  0)  +  0(e2)  ,  ( 

xy 

/,n  da(0)(x,  0) 

°ns(x‘  X)  “  axy)(x’  0)  +  °xj)(x’  0)  "  X(x)  “ ^ -  + 

+  X«(x)[a^)(x,  0)  -  «£>(*,  0)]  +  0(e2)  ,  ( 

yy  xx 


un(x,  X) 


uy°\x,  0)  +  uy^(x,  0)  +  X(x)syy^(x,  0) 


-  X'(x)u<0)(x,  0)  +  0(e2)  . 


Finally,  substituting  eqns.  (9)-(ll)  and  ( 15)— ( 16)  into  the  boundary  value 


problem  formulated  in  Section  3.1  (eqns.  (l)-(4)>,  taking  into  account  (5)-(8) 
and  separating  zero  and  first  order  terms,  we  find  that  and 


should  be  the  solution  to  the  following  boundary  value  problem 
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3.  FORMULAE  FOR  THE  STRESS  INTENSITY  FACTORS 

Following  Cotterell  and  Rice  [7],  let  u>  be  the  angle  of  the  crack  tip 
at  x  =  a  ,  given  by  u  =  \'(a)  to  first  order.  The  normal  (o^)  and  shear 
(oraJ)  stresses  acting  along  the  prolongation  of  the  crack  at  a  small  dls- 
tance  r  from  the  tip  at  x  =  a  are  obtained  by  setting  \  =  wr  +  ©(e’)  = 
u(x  -  a)  +  O(e^)  into  equations  (15)  and  (16).  So, 


[«r« 


S 


<J  -  a  (x,0)-t»>(x-a)  - Hr — 

aw  yy  3x 


-2<ja  (x,0)+o'  (x,0)+0(e  )  , 

yy 


o) 


a  ■  a^^(x,0)Hi)(x-a)  — - +ii>[o^(x,0)-o^(x,0)]+<4^(x,0)+0(e2)  . 

ru  xy  *  3x  yy  ’  xx  ’  xy 


Then,  the  stress  intensity  factors  can  be  calculated  as 


Kt  -  lim.  </2wr  a  )  =  4°^  +  K*1*  +  4^  +  0(e2)  >  (2 

I  rt+  axo  i  Iu  l  ’ 

r-*-0 

K-  -  &  ^  •«>  -  Kn  ♦  4\L  ♦  4?  *  0(c2)  , 

where  *  4^  an<*  4l^  are  t*ie  stress  intensity  factors  for 

the  zeroth  (eqns.  (5)-(8))  and  first  order  (eons.  (18)-(22))  problems,  and 


=  -w/H  lim+[(x  -  a)3/2  g-A-! - -  +  2(x  -  a)1/2  o^}(x,  0)] 

x->a  y 


=*  -u/2x  lim  {(x-a)3/2  ^  +  (x-a)1/2[o(0) (x,0)  -  o(0)(x,0)]}  . 

IIw  .  +l  9x  Lyy  xx  J) 


Using  the  last  two  equations  and  a  Williams  [11]  expansion  for  the  near 
crack  tip  stress  field,  we  can  show  that 


_  1  (0) 

2  KII  ’ 


and 


i  .  K<°> 


From  the  formulation  of  the  first  order  problem  (eqns.  (18)-(22)),  it  is 


clear  that  this  can  be  considered  as  the  superposition  of  the  following  two 
problems;  problem  (i)  with  a  prescribed  normal  dlsplacment  and  zero  shear 
traction  on  the  sliding  portions  of  the  crack  and  with  the  rest  of  the  crack 
traction  free,  and  problem  (ii)  with  a  prescribed  shear  traction  and  zero 
normal  traction  everywhere  on  the  crack  face.  Thus,  knowing  the  solution  of 
the  zeroth  order  problem  (eqns.  < 5)— (8))  and  having  determined  the  sliding  and 
open  portions  of  the  crack,  we  can  determine  from  the  solution  of 

problem  (1)  mentioned  above. 

As  far  as  is  concerned,  it  is  obvious  that  only  the  prescribed 

shear  tractions  at  the  crack  surface  of  problem  (ii)  mentioned  above  that  have 
opposite  directions  on  the  upper  and  lower  surfaces  of  the  crack  have  a  non¬ 
zero  contribution  to  .  With  the  definition 

A(x)  -  j  [A(x,  0+)  +  A(x,  0")]  , 


is  known  (e.g.,  [12])  to  be 


1 

Zita 


jV‘>(,) 

j  Xy 
-a 


/ 


a  +  x 
a  -  x 


dx  , 


(27) 


where,  according  to  (21), 


xy 


(x) 


^°(x)  +  \(X)  4-  r^0)(x)  -  i^0)<x)i  - 

yy  dx  -  xx  '  xy  ' 

-  X' (x)[  (1  +  2*i2)^<°)(x)  -^(x)]  . 

yy  xx 


(28) 


On  the  other  hand,  it  is  possible  that  the  stress  field  of  the  zeroth  order 
problem,  ,  has  the  characteristic  —  elastic  singularity  at  several 


points  in  the  interval  |x|  <  a  ;  since  derivatives  of  a  with  respect 

to  x  are  involved  in  the  formula  for  (eqn.  (28)),  non-integrable 

xy 

singularities  will  appear  in  eqn.  (27).  To  overcome  this  difficulty,  we 
assume,  for  the  moment,  that  the  stress  components  cr^^(x,  0)  are  all 


bounded  and  differentiable  with  respect  to  x  in  the  interval  ]x|  <  a  ;  this 

makes  c/ ^ (x,  0)  also  bounded  on  the  crack  face.  In  the  case  where 
xy 

o£^(x,  0)  are  singular  at  some  point  in  the  interval  |x|  <  a  ,  the  singu¬ 
larities  are  removed  by  replacing  a^?\x,  0)  by  bounded  functions  that 
reduce  continuously  to  zero  (or  any  other  value  that  makes  o^^(x,  0) 
continuous)  over  distances  closer  than  a  small  distance  6  to  the  point  where 
the  singularities  appear.  Later  it  is  shown  that  it  is  possible  to  let  6 
tend  to  zero,  i.e.,  effectively  to  remove  the  restriction  of  bounded  and 
differentiable  o^j\x,  0)  . 

We  return  now  to  the  calculation  of  K  •  With  the  above  continuity 

assumptions  on  we  can  integrate  by  parts  eqn.  (28)  to  find 


k(t;}  -  -  jl  /  °  -  x '[ ( i  +  2pz)^o)  -  ^0)]  + 

11  /s  -a  yy  yy  504 


+  [-*'  ♦{x,<.>irc£> -,«<»)»  /fii  ax 


-  —  /8  (^0)  -  ^0)  H  7  X’(a)  -  -a-?L.-t.(a  -  x),x:,(a),  /  a_LJ£  dx 
, —  •''■xx  'xy  '  -  2  ,  n2  a  +  x 

/n a  -a  (a  -  x) 


It  should  be  noted  that  aX(x)  +  (a-x)X’(a)  =  -r—  (aX(x)  +  (a-x)X’(a)]  =*  0 

dx 

at  x  =  a  ,  so  there  is  no  divergence  at  the  upper  limit  of  the  second  inte¬ 
gral  in  (29).  It  can  also  be  seen  that  an  integrable  singularity  can  exist  in 
0^j\x)  ,  provided  it  is  not  at  x  =  ±a  ,  as  was  also  noted  in  (7).  Specifi¬ 
cally,  in  terms  of  our  earlier  discussion,  6  can  be  shrunk,  to  zero  and  in 


11 


that  limit,  the  result  of  equation  (29)  for  Kjj.  approaches  the  result 
obtained  by  Inserting  directly  into  (29)  the  singular,  actual  .  Such 

considerations,  based  essentially  on  the  fact  that  the  final  result  of  equa¬ 
tion  (29)  for  contains  ^^(x)  only  (and  not  derivatives  with  respect 

—(0) 

to  x),  allow  us  to  conclude  that  (29)  is  valid  for  all  integrable  (x) 

(i.e.,  not  necessarily  bounded  or  continuous). 

We  mention  again  that  part  of  the  solution  to  our  problem  is  the  finding 
of  the  sliding  portions  of  the  crack.  Having  found  the  sliding  portions  of 
the  crack  and  the  solution  of  the  zeroth  order  problem,  we  can  proceed  to 
solve  the  first  order  problem  and  use  the  formulae  given  in  this  section  to 


find  the  first  order  correction  to  the  stress  intensity  factors. 


4.  THE  PROBLEM  OF  THE  KINKED  CRACK 

A  particular  case  of  the  curved  crack  is  the  kinked  crack  shown  in  Fig. 
2.  The  shape  of  the  kinked  crack  is  given  by 


r  mb 
b  -  a 

(x  +  a) 

for  -a  <  x  <  -b  , 

X(x)  =  mx 

for  | x  j  <  b  , 

mb 

(x  -  a) 

for  b  <  x  <  a  . 

b  -  a 

In  this  case,  w  =  \'(a)  +  O(e^)  =  t~ —  +  O(e^)  . 

D  i 

Following  our  previous  discussion,  we  assume  that  both  m  and 
are  0(e)  ,  which  is  equivalent  to  assuming  that  X’(x)  is  0(e)  .  We 
mention  agaLn  that  we  are  concerned  with  the  case  where  both  Oy“  and  oXy 


are  negative  and  so,  it  can  be  assumed  that  the  portion  of  the  crack  in  the 
interval  |x|  <  b  remains  closed  during  the  application  of  the  load.  Thus, 
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the  sliding  portion  of  the  crack  is  the  interval  |x|  <  b  and  the  open  por¬ 
tions  are  the  intervals  b  <  |x|  <  a  .  Of  course,  there  is  a  possibility  that 
the  applied  stresses  o^j  and  the  orientation  of  the  kinked  crack  are  such 
that  the  whole  crack  remains  closed  and  does  not  slide.  To  check  whether  this 
happens,  we  can  solve  the  problem  assuming  that  the  crack  opens  in  the  inter¬ 
val  b  <  |x|  <  a  ;  if  the  calculated  at  the  tips  of  the  kinks  is 

negative,  the  tip  of  the  kinks  remain  closed  and  the  assumption  that  the  crack 
opens  in  the  whole  intervals  b  <  |x|  <  a  is  in  error. 


4.1  Solution  of  the  zeroth  order  problem 

The  zeroth  order  problem  can  be  considered  as  the  superposition  of  the 
four  problems  shown  in  Fig.  3,  where  F(x)  is  the  distribution  of  the 
0y^(x,  0)  stress  component  of  problem  no.  1.  We  note  that  for  problem  no.  4 
the  shear  stress  on  the  crack  face,  cf^^(x,  0)  =  pF(x)  ,  opposes  the  relative 

xy 

sliding  of  the  crack  faces.  The  quantities  of  interest  for  each  of  the  four 
problems  mentioned  above  are  given  in  the  following.  In  the  solutions  pre¬ 
sented  in  the  rest  of  this  section,  conditions  of  plane  strain  are  assumed;  in 
order  to  get  the  plane  stress  solutions  we  simply  replace  v  by 


1  +  v  * 


4.1.1  Problem  no.  1 


The  solution  to  tnis  problem  has  been  given  by  Erdogan  [5]  and  is  as 
follows 


a(0)(x,  0)  = 

yy 


yy 


2  E(k)  2 


//,  2  2W  2  2. 

/(b  -  x  )(a  -x) 


K(k) 


-  x 


]  =  F(x)  , 


o<°> 

XX 


(x, 


0) 


F(x) 


9 


a <°> 

xy 


(x. 


0)  =  0 


I 


(30) 


u^"(x,  O')  «.  u^"(x,  0  )  , 


€<0)(X,  0+)  -  e(0)(x,  0  ) 

yy  yy 


for  x  <  b  ; 


c/^(x,  0)  »  -o  °° 
xx  yy 


for  b  <  x  <  a  , 


where  K(k)  and  E(k)  are  Che  complete  elliptic  integrals  of  the  first  and 

2  2 

second  kind  respectively  and  k=*/l-b/a  . 

Also,  the  stress  intensity  factors  for  this  problem  are 


K(0)  =  V/lta  f!  -Ml 

KI  k  K(k)  J  * 


and  -  0  . 


Problem  no.  2  consists  of  a  plane  strain  tension.  The  solution  to  this 
problem  is  quite  obvious;  therefore  we  proceed  to  problem  no.  3. 

4.1.2  Problem  no .  3 

The  solution  to  this  problem  is  known  (e.g.,  [12])  to  be 


(0),  «°  2x 

i  (x,  0  )  =  +  a  — - -  , 

w  '  vw  - — —  * 


,  2  2 
/a  -  x 


c(0)(x,  0*)  =  a(0)(x,  0*)  =  0  , 
yy  xy  *  ’ 


e  (x,  0‘)  -  ±v  - - - 

yy  G  ,  2 — 2 

/a  -  x 


and  u^(x,  0*)  =  ±(1  -  v)  -jjp-  ■/&"  -  x^  for  |x|  <  a  , 


where  6  Is  the  shear  modulus  of  the  material. 


In  addition. 


l"’  -  0  . 


and  =  a  ™  /ica 

II  xy 


A. 1.3  Problem  no.  4 

The  solution  to  this  problem  is  derived  In  Appendix  1  and  is  as  follows 


o<0>(x,  0*)  -  "  2* 

XX  ’  yy  -5 - r 

/a  -  x^ 


°yy)(x’  0)  "  0  ’ 


cr^Cx,  0)  =  pF(x)  , 
xy 


e<0)(x  0*)  -  +v  —2 S2. _ 5 _ 

Eyy  °  ^  G  - - x  ’ 

/a  -  x 


u^(x,  0*)  =  +(1  -  v)  — /a^  -  x^  for  Ixl  <  b  ; 
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o(0)(x,  0*) 
xx 


rx>  , 

+2pa  I 


2  2  E(k) 


2  2  2 
b  l(a^  -  x  ) 


,  2  2 
/a  -  x 


for  b  <  x  <  a  . 


Also, 


K^0)  -  0  , 


c<0>  -  ft  _  k  _EOO  i 

II  11  K  K(k)  J 


4.1.4  Superposition 


Superimposing  the  solutions  of  the  four  problems  shown  In  Fig.  3,  we  find 
the  quantities  of  interest  of  the  solution  to  the  zeroth  order  problem  to  be 


— ( 0 )  00  00 

o'  '(x)  -  F(x)  -  o  +o  , 
xx  yy  xx 


^0)  -  F(x)  , 

yy 


^(x)  =  pF(x)  , 
xy 


E(0)(X,  0+)  -  e(0)(x,  0*)  =  2V  - * - 

yy  yy  g  ^ ^ 

/a  -  x 


u(0)(x,  0+, 

X 


O  -  \LO  -r - r* 

'  (x,  0-)  -  2(1  -  V)  -g— — TC  U1  -  x2 


for  x  <  b  ;  and 


- (0)  x  x  00  °» 

o'  '(x)  =  -a  +  a  , 
xx  yy  xx 


(36) 
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^0)(x)  -  ^<0)(x)  -  0 


for  b  <  jx  <  a  . 


Also,  the  stress  intensity  factors  for  the  zeroth  order  problem  are  given 


r(°)  _  gyyy/rc*  r,  _  E(k)  -I 

h  k  Li  K(k)  J  * 


4V  -  -  O  Aa  +  ^0) 


4.2  Stress  Intensity  factors  for  the  first  order  problem 

As  discussed  In  Section  3,  the  first  order  problem  can  be  considered  to 
be  the  superposition  of  the  two  problems  shown  in  Fig.  4. 

The  mode  I  stress  intensity  factor  is  determined  by  solving 

problem  no.  1  in  Fig.  4,  which  is  actually  the  problem  of  the  opening  of  a 
finite  crack  by  a  rigid  wedge.  The  general  solution  to  this  problem  has  been 
given  by  Markuzon  [14].  Taking  into  account  eqn.  (22)  and  the  solution  of  the 
zeroth  order  problem  derived  in  the  previous  section,  we  find  that,  for  our 
particular  case,  the  shape  of  the  wedge,  h(x)  ,  is  given  by 


<»  <»  2  2 

2h(x>  -  u(1>(x,0+)-u(1)(x,0-)  -  2m  JSL^EL  Qz<0jL.1L.  ,  1*1  <  b  .  (40) 


/„2  -  x2 


Using  the  above  formula  for  the  shape  of  the  wedge  and  Markuzon's  [14] 
solution,  we  find  the  mode  1  stress  intensity  factor  to  be  (see  Appendix  2) 


.,(1)  ,  ®  “x  / —  m  k  -  (2  -  v)k  +  (l  -  2v  +  2C)k  +  v 

Kr  '  -  -(c  -  pa  )  vna  - - - - - - - -  , 

i  *y  yy  l  -  v  2l  2  ’ 


(41) 


•K(k)k  +  [(1  +  2v)K(k)  -  2E(k)]k  +  2v 


We  proceed  now  to  the  calculation  of  the  mode  II  stress  intensity  factor 

for  the  first  order  problem,  >  which  can  be  determined  either  by  solving 

problem  no.  2  in  Fig.  4  or,  equivalently,  using  eqn.  (29).  We  note  that 

cr^^(x)  ,  o^^(x)  and  o^^(x)  all  have  the  characteristic  —  singularity 

xx  yy  xy  . 

(iy  r 

at  x  =  ±b  ;  but  as  discussed  in  Section  3,  eqn.  (29)  for  Kjj.  can  still  be 
used,  provided  that  the  singularities  are  integrable,  which  is  indeed  the 
case.  So,  using  the  solution  of  the  zeroth  order  problem  derived  in  the  pre¬ 
vious  section  and  applying  eqn.  (29),  after  some  lengthy,  but  straightforward, 
integrations  we  find 

„  ( 1 )  „(1)  .  /  •  00 -v  / —  mb  t .  2a  .  — 1  b  \  2  <*>  . — 

Kt;  *  pKT  +  (cr  -a  )  /xa  -r -  [1 - r-  sin  —  ]  +  2li  mo  /xa  + 

II  ^  I  yy  xx  b  -  a  v  xb  a  1  K  yy 


where  is  given  in  eqn.  (40). 


5.  DISCUSSION  AND  COMPARISON  WITH  THE  EXACT  SOLUTION 

The  obtained  asymptotic  solution  for  the  stress  intensity  factors  at  the 
tips  of  a  kinked  crack  is  of  the  form 


K(°>  3  K<0>  +  K(l>  +0.  2.  _ 


K 
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KII  "KI?)  +J“  KI°)  +  ^  +  0(e2)  » 

where  and  Kj^  are  given  in  equations  (38),  (39), 

(41)  and  (42)  respectively.  Equations  (41)  and  (42)  show  that  and 

depend  on  Poisson's  ratio  v  .  This  means  that  they  have  different 
values  under  plane  strain  or  plane  stress  conditions.  On  the  other  hand,  it 
is  known  that,  in  the  absence  of  body  forces,  the  stress  intensity  factors  for 
traction  boundary  value  problems  are  the  same  under  plane  strain  or  plane 
stress  conditions  and  independent  of  the  elastic  constants.  The  reason  that 
v  enters  the  expressions  of  the  asymptotic  solution  for  the  stress  intensity 
factors  is  that  the  displacement  field  of  the  zeroth  order  problem,  which 
depends  on  v  ,  was  used  to  formulate  the  first  order  problem.  As  a  result, 
the  in-plane  components  of  the  first  order  correction  to  the  stress  field  and 
the  corresponding  corrections  to  the  stress  intensity  factocs  depend  on  v  . 
However,  numerical  calculations  of  and  show  that  their  depen¬ 

dence  on  v  is  very  weak  and  that  their  values  for  plane  strain  and  plane 
stress  conditions  are  practically  indistinguishable,  which  validates  the  first 
order  correction. 

Next,  we  apply  our  results  to  the  problem  of  a  infinite  plate  containing 
a  kinked  crack  oriented  at  36°  to  the  overall  compression  (Fig.  5).  The  exact 
solution,  given  in  Ref.  [3],  and  the  asymptotic  results  for  Kj  are  plotted 
in  Fig.  6  versus  the  angle  between  the  straight  crack  and  its  out-of-plane 
kinks,  9  ,  for  several  values  of  the  ratio  of  the  length  of  the  kink,  X  ,  to 
the  length  of  the  straight  crack,  c  .  In  general,  the  region  of  accuracy  of 
the  asymptotic  solution  depends  on  both  X/c  and  9  ,  because  the  values  of 
m  and  co  in  our  analysis,  which  must  be  small  for  the  asymptotic  solution  to 


be  valid,  depend  on  both  1/c  and  0  .  But,  roughly  speaking,  the  asymptotic 
solution  is  seen  to  be  accurate  for  values  of  9  up  to  about  20°.  Unfortu¬ 
nately,  the  values  of  Kjj  for  small  values  of  6  are  not  given  in  Ref.  [3]; 
so,  comparisons  of  the  asymptotic  result  for  Kjj  with  the  exact  solution 
were  not  possible. 

6.  CLOSURE 

A  first  order  solution  has  been  obtained  for  the  stress  intensity  factors 
at  the  tips  of  the  kinked  extension  of  a  sliding  crack.  The  validity  of  the 
asymptotic  solution  is  limited  to  kinked  cracks  with  small  deviations  from 
straightness.  There  are  several  situations  where  this  deviation  is  indeed 
small.  As  an  example,  consider  the  case  of  a  glass  plate  or  a  rock  block  con¬ 
taining  several  small  cracks  at  different  orientations.  Under  the  application 
of  a  compressive  load,  the  cracks  with  an  angle  to  the  direction  of  compres¬ 
sion,  y  ,  greater  than  yc  =  tan  1  -jj  will  remain  closed  and  only  those 
with  y  <  yc  can,  possibly,  slide  and  propagate.  It  is  also  known  ([l]-[4]) 
that  these  cracks  tend  to  propagate  towards  the  direction  of  compression.  So, 
if  the  coefficient  of  friction,  p  ,  is  very  high  (which  makes  yc  small)  the 
crack  propagation  will  create  kinked  cracks  with  small  deviations  from 
straightness.  For  situations  like  these,  the  asymptotic  results  can  be  used 
to  determine  the  stress  Intensity  factors  and  to  make  predictions  for  the 
direction  of  further  propagation.  In  addition,  fatigue  due  to  non¬ 
proportional  loads  can  cause  the  development  of  cracks  that  are  not  straight 
and  are  partially  closed,  although  open  at  the  tip.  For  cases  where  the 
deviation  from  the  straight  line  is  small,  the  methods  devised  here  can  be 
used,  although  a  criterion  for  determining  where  the  closed  portions  lie  would 


have  to  be  developed. 
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APPENDIX  1 

The  problem  no.  4  In  Fig.  3  Is  formulated  in  terms  of  the  complex 
potentials  $  and  4*  of  Muskhellshvili  [10).  The  stresses  and  displacements 
can  be  expressed  as 

°xx  +  °yy  "  2l*’(z>  +  »  (43) 

Oyy  -  a ^  +  210^  =»  2[z<J>"(z)  +  4>’(z)]  ,  (44) 

2G(ux  +  itiy)  -  <$(z)  -  zifr'U)  -  4>(  z)  ,  (45) 

where  k  ■  3  -  4v  for  plane  strain  and  <  =*  (3  -  v)/(l  +  v)  for  plane 
stress,  the  overbar  denotes  the  complex  conjugate  and  prime  stands  for 
differentiation  with  respect  to  z  *  x  +  ly 
Introducing  the  analytic  function 

fl(z)  <=  z<b'(z)  +  <i»(z) 

eqns.  (43)-(45)  can  be  written  as 

°xx  +  °yy  “  2l*'(z)  +  <5>'(z)]  »  (46) 

a  -  a  +  2ic 3  =  2[(z  -  z  )<*>"(  z )  -  Q(z)  -  $'(z)]  ,  (47) 

JJ  *7 

2G(u  +  iu  )  -  K<J>(z)  -  (z  -  z)  <t>'(z)  -  Q(z)  . 

x  y 


(48) 


For  a  prescribed  shear  traction  oxy(x,  0)  along  the  crack  face,  it  is 
known  (e.g.,  [12])  that 


♦  ’(z)  -  Q'(z) 


— -  f  a  (x,  0)  /a2  - 

- -  -  J  vu  ' 


o/2  2  -a 

2x  /z  -  a 


2  2  dx 

_  x  - 

x  -  z 


In  our  problem 


axy(x,  0) 


pF(x)  for  | x |  <  b  , 

0  for  b  <  |x|  <  a  , 


where  F(x)  is  defined  in  eqn.  (30).  Substituting  (49)  into  (50)  and 
carrying  out  the  integration  we  find 


<fr’(z)  -  Q’(z)  -  -  y  poy”{[z2  -  a2 


✓(z2  -  a2)(z2  -  b2)  /z2  -  a2 


Finally,  using  eqns.  (46)-(48)  and  the  definition 


KI  +  1KII  “  lim+  ~  a)  [<Jyy(x»  °>  +  i0xy^x’ 


we  find  the  results  shown  in  Section  4.1.3. 


APPENDIX  2 


The  solution  to  the  problem  of  the  opening  of  a  finite  crack  by  a  rigid 
wedge  has  been  given  by  Markuzon  [14].  In  terras  of  Muskhelishuill's  [10] 
complex  potentials,  the  solution  is  shown  to  be 


»  « 


m 


$'(z)  =•  Q'(z) 


2G  JL _  fb  dh  ,  .  dx  Co 

x(<  +  1)  X(z)  !  dx  X  x  -  z  X(z)  * 


where  h(x)  is  the  function  determining  the  shape  of  the  wedge  of  length  2b 

2  2  2  2 

(see  Fig.  4),  X(z)  =*  /(a  -  z  )(b  -  z  )  and  the  constant  CQ  Is  determined 


from  the  equation 


1  / 


r  t  dh  2  2x/,2  2»  dt  -i  , 

f  [  3t  /(a  ~  c  )(b  '  c  >  — IT ldx  ‘ 


b  j.  2  2W  2  ,2.  -b  A 

/(a  -  x  )(x  -  b  ) 

k  +  1  K(k)  „ 

~~2G  ~a  Co  -  _h(b)  * 

As  discussed  in  Section  4.2,  the  shape  of  the  wedge  for  our  problem  is 


given  by 


o  *  -  pa  "  ,  v  2  2 

1.00  ry  a  -  »)»  .  U  <  b  . 


J  2  2 

/a  -  x 


Substituting  (53)  into  (52)  and  carrying  out  the  integrations,  we  find 


4>'(z)  =  Q'(z)  =  - 


2(1  -  v) 


CO  CO 

0  -  po 

xy  yy 


„  2  2W  2  ,2. 

/(z  -  a  )(z  -  b  ) 


z/z2  -  b2  -  a/a2  -  b2  .  z2  -  z/z2  -  b2  2va2  +  b2 


K(k)k  +  [(1  +  2v)K(k)  -  2E(k)]k  +  2v 
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Finally,  using  eqns.  (46),  (47)  and  (51)  we  find  the  stress  intensity  factor 
to  be 


lict  )  /ita 

yy 


k3  -  (2  -  v)k2  +  (1  -  2v  +  2C)k  +  v 


1  -  v 


Fig.  1.  Infinite  plate  containing  a  curved  crack. 

Fig.  2.  Infinite  plate  containing  a  kinked  crack. 

Fig.  3.  Superposition  used  in  the  solution  of  the  zeroth  order  problem. 

Fig.  4.  Superposition  used  in  the  solution  of  the  first  order  problem. 

Fig.  5.  Infinite  plate  containing  a  kinked  crack  oriented  at  36°  to  the 
overall  compression. 

Fig.  6.  Stress  intensity  factor  at  the  tips  of  the  kinked  crack  shown  in  Fig 


■  Exact  solution 
Asymptotic  solution 


